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A quantum spin model representing tautomeric mutation is proposed for any DNA 
molecule. Based on this model, the quantum mechanical calculations for mutational rate 
and complementarity restoring repair rate in the replication processes are carried out. A 
possible application to a real biological system is discussed. 



1 



1 Introduction 



It is well known by now that in the cellular environment, DNA molecules are not absolutely 
stable and invariant. Rather each base pair in a DNA double helix has a certain probability of 
undergoing change. Such events occurring within a gene are known as "gene mutations" [1] . 
Moreover DNA can be regarded as a system undergoing a dynamic competition between the 
chemical processes leading to new mutations and the cellular repairing processes correcting 
these premutational events [1,2]. 

Here we would like to present a possible model and calculational scheme corresponding 
to a non-enzymatic partial repairing process leading to "complementarity restoration" for 
a special kind of point mutation known as "tautomeric mutation" [1,2]. The origin of this 
type of mutation is quantum fluctuation occurring spontaneously or at times initiated by 
the presence of mutagens (base analogs). In this mutational process each nucleotide (viz. 
A,T,C,G) can fluctuate between the two states viz. "keto" and "enol" driven by "proton 
tunnelling" [2]. Tautomerization is known to affect approximately 0.01 percent of the DNA 
bases (nucleotides) [2]. Out of these two states the most commonly occurring keto states 
obey the conventional Watson-Crick base pairing rule (BPR), whereas the rare enol states can 
have various anomalous pairings [1,2]. During replication the presence of this enol form (if 
present in any of the strands) leads to anomalous base pairing and hence causes "tautomeric 
mutation" (TM). After the completion of more than one round of replication however, the 
original base pairs, involved in normal pairing on the two strands, may change over to a 
new normal base pair leading to "base pair transformation" (BPT); Or else they may come 
back to the original unchanged base pair itself. This phenomenon is called "complementarity 
restoration" (CR). 

In this communication we would like to calculate the CR rate in such a TM process. 
Our calculation is quantum mechanical and based on a quantum spin modelling of the DNA 
configurations. The model we are adopting is an analog of a 2-legged Ising spin ladder [3-6] 
involving two kinds of "magnetic atoms" in the presence of a "fictitious transverse field", 
with only inter-strand interactions. This inter-strand antiferromagnetic interaction of the 
Ising type operating between the "spins" on the "like atoms" only, tries to preserve the 
"Watson-Crick base pairing rule (BPR)" or the "complementarity" on the two strands of 
the DNA. Furthermore, these two kinds of atoms themselves are also represented by another 
kind of "pseudo-spin" | objects experiencing only an inter-strand "pseudo-ferromagnetic 
interaction" of Ising type, to ensure the presence of like atoms on the two strands. The 
intra-strand process modelled on the assumed coupling of both "spin" and "pseudo-spin" 
degrees of freedom to a fictitious transverse field, tends to cause "mutational fluctuations" 
leading to violation of the BPR. We set up the equations for calculating the rates for both 
BPR violating point mutations as well as CR processes during the replication, for any ar- 
bitrary DNA sequence containing N nucleotides on each strand. From this we estimate the 
efficiency of the "secondary repair mechanism" for a typical DNA molecule with some known 
parameters. We suggest a possible system and experiments to test our theoretical prediction. 

It may be mentionworthy that simple transverse Ising model has been used to study 
proton tunnelling in a double well potential configuration by many workers [3,4]. Our model 
and approach may therefore be looked upon as a non-trivial extension of this simple idea to 
the biological processes at the molecular level. 
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2 Mathematical Formulation 



We model "point mutations" of tautomeric type by a quantum spin model. As mentioned 
in the last section, the model we employ is analogous to a transverse field Ising model of a 
2-legged ladder [3-7]. More explicitly, it is given as:- 

nj n/inter—strandrung . nj intra— strand /i\ 

n — n Ising ~r n Trans \ l ) 

Where, 



Rising — ^2 ^i 1 ' 2 ^ ' i,M,S]l a * ' i,N,s';2&MN&s,-s' (2) 
i,M,N,s,s' 

with the variables M and N both representing varieties of atoms X and Y and a z taking 
values s and s', the "spin values" corresponding to these atoms. Both s and s' can assume 
numerical values | and The indices 1 and 2, occurring as superscripts and subscripts 
in A and a respectively, represent the two strands. The inter-strand coupling constant Aj 1 ' 2 
connecting the nucleotides on the ith pair of sites on the two strands, is assumed to be -ve. 
The nucleotides A,T,C and G can then be identified with the states (X, |), (X, ^-), (Y, |) 
and (F, respectively. Moreover, X and Y are themselves represented as pseudo-spin "up 
state" and pseudo-spin "downstate" respectively. 

We have assumed in the above model that the intrinsic coupling is purely inter-strand and 
that the intra-strand intrinsic coupling is zero. This inter-strand coupling is antiferromag- 
netic for the spins but ferromagnetic for the pseudo-spins, as is obvious from the structure 
of the above hamiltonian. This tries to preserve the BPR allowed nucleotide pairings on 
the two strands. Each strand is individually and internally "paramagnetic" (uncorrelated in 
both spin and pseudo-spin) but they are mutually both "spin and pseudo-spin correlated 
(paired)". In this model the BPR violating mutations are generated by the transverse field 
(assumed to be external) acting on each strand. This part of the Hamiltonian is given by:- 

1 1 

yintrastrand = ^ -H\[(a^ + <J^) sp XI psp + I sp X(af l + <7 t ~,)psp + -(of, + (7~;) sp X((J^ + CT t ~,)pap] 
id 

(3) 

where sp stands for "spin" and psp stands for "pseudo-spin". The cr's are the Pauli matrices 
and / is an unit matrix of 2X2 dimension. It should be highlighted that our model differs 
from the conventional transverse Ising model in two major aspects:- (i) the transverse field 
H\ couples to both spin and pseudo-spin locally and (ii) both the couplings may be present 
simultaneously. Moreover the site dependence of the transverse field ensures the different 
amplitudes for various mutational channels. 

Let us assume that the magnitude of Aj is much greater than that of Hi. We then attempt 
to treat the 2nd term in the above Hamiltonian, as a perturbation and use the ground state of 
the 1st term (in the Hamiltonian) as the reference state. The ground state is a configuration 
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which obeys BPR. We confine ourselves to the case of zero temperature in the absence of any 
external stimuli (like UV radiation). In our scheme the order of perturbation corresponds 
to the order of replication of the DNA. The perturbation calculation will have to be based 
on the degenerate perturbation theory, as the ground state has a high degree of degeneracy, 
viz. 4^, where iV is the total number of nucleotides in any one strand of the DNA molecule. 
Let l^j) and \ipj) be the two degenerate ground states and \if) ex ) be the excited state. Thus 
for our double-stranded DNA molecule, is a "parent" molecule state, \ip ex ) represents 
one of the "daughter" molecule product states and \ipj) is the "grand daughter" molecule 
state. The daughter state is obtained from the parent state after the 1st round of replication 
causing some mutation. Again the grand daughter state is obtained from the parent state 
after the 2nd round of replication and which is one of the doubly mutated states degenerate 
with the parent state. Thus after even number of replications (even orders of perturbation), 
the successive mutations can sometimes be mutually compensating to get back a new state 
satisfying BPR but related to the original parent state by a BPT. The other possibility is to 
come back to the original state itself. The odd number of replications however, always gen- 
erate some BPR violating defective configurations. Besides for a real DNA molecule, even 
after the 1st round of replication there are normal (unmutated) configurations produced 
which will also continue to lead to the original parental state after the next (2nd) round of 
replication as well. Our model however does not include these normal replication processes. 
We only model the processes which involve both mutation and replication. Moreover we only 
allow single (point) mutations one at a time and neglect the doubly (and higher) excited 
states in our calculations. 

We would like to study here the lowest order processes viz. the second order replication for 
calculating the "complementarity restoration" rate (CRR) and the first order replication for 
calculating the tautomeric mutation rate (TMR). For these calculations, we will follow the 
standard procedure from quantum mechanics [8]. 

We would need the following two quantities for our calculations:- 

Fl(z,ex) = mn I T 3 rans \r x ) (4) 

and 

F2(i,j) = £ (WEanslr^i^lKljilj) (5) 
ex 

where is the "parent state", \if) ex ) is a "daughter product state" and \tpj) is a "grand 
daughter product state" , as have been discussed earlier. The symbols i, j and ex denote the 
indices corresponding to the two degenerate ground states viz. |^) and respectively and 
the excited state \ip ex ). Both \ipi) and belong to the manifold containing 4^ degenerate 
ground states, where iV is the total number of nucleotides on any one strand of the DNA 
molecule. Again for a chosen ground state, the degree of degeneracy of the lowest excited 
states (containing single mutation) is 2N, as can be seen by combining mutation with repli- 
cation. Moreover for a definite \ipi) the grand daughter state \ipj) has only 2 choices for the 
complementarity restoration processes, corresponding to each mutated daughter state \ip ex ). 
As discussed earlier, one of these two states corresponds to the original parent state and the 
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other one a state obtained through BPT. This leads to the total number of distinct grand 
daughter states as 2N + I. 



These quantities F\ and F2 will be the important ones in the calculations involving the 
first order and the second order processes respectively. The rate for the base pair transforma- 
tion can be calculated by making use of the Fermi Golden Rule [8], as the parent state and 
the grand daughter state are degenerate. The calculation for BPR violating mutational rate 
is however more complicated, as the process does not conserve energy strictly. The deviation 
from the energy conservation is however times the absolute magnitude of the ground state 
energy. Thus for a macroscopic size of the DNA molecule, even the BPR violating single 
mutation conserves energy approximately. Therefore we can estimate the rate of this process 
by applying Fermi Golden Rule once again. We derive an approximate expression for the 
rate in this case. 

-iTotai [2nfY,ij pi\Fl{iJ = ex)\ 2 5(uj fi ) 

T NBPR — ^ W 

Here, pi is the probability of occupation of the initial (chosen) ground state with index i, 
Ufi is simply / is the index of the final (lowest excited) state and AE is the energy 

difference between the excited state and the ground state. As has been pointed out before, 
for the single mutation F AE tends to zero. It can be seen very easily that in general pi is 



ground 



independent of i and is equal to 4 N . Similarly, the rate for CR is given by 

_ lT otai_ {2ir) 2 pEij\F2{i,j)\ 2 8{u ij ) 



BPR 



h 



(7) 



From the above two expressions, we see that the rates for both the processes contain the 
matrix element of Ti.^ rans between a ground state and an excited state configuration. Thus 
it should be possible to relate the two rates and thereby estimate the repairing efficiency 
(RE). By taking the ratio of the CR rate and the TM rate, with the proper degeneracy 
factors, we get the expression for RE. Assuming an uniform external field H for simplicity, 
the expression for RE becomes equal to 2H 2 . 



3 Calculations and Results 

As an example we now estimate the magnitude of the fictitious transverse field H, the 
simulator of mutation introduced in our model, from the available data for spontaneous 
point mutation corresponding to a typical DNA molecule. It is found that in the absence of 
mutagenic treatment the effective spontaneous mutation rate (after primary repair processes 
or "proof reading") is about 2X10" 11 per second per nucleotide for a well studied organism 
(a fungus) named "Neurospora" [1]. Since in this case the parental DNA configuration and 
the mutation channels are all fixed, assuming that all the mutations referred above are of 
tautomeric type having taken place in the 1st round of replication itself, we get the following 
equation involving the rate:- 

H 2 n = NX2X1Q- 11 (8) 
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where n is the number of mutation channels (mutated daughter states) allowed after the 1st 
round of replication, for this species. As mentioned before, n is proportional to N and in fact 
is equal to 2N. Thus the magnitude of H becomes approximately 3X10 -5 in the suitable unit 
and is essentially independent of N. This leads to the magnitude of the secondary repairing 
efficiency of the order of 10~ n . To be more precise, according to our calculation about 2 out 
of every 10 11 point mutations are compensated (complementarity is restored) through the 
secondary repairing processes arising from quantum fluctuation itself, for Neurospora. Thus 
this number will change with different choices of the species only if the rate of mutation 
changes [9]. 

The results will be very different however, if we start from an arbitrary DNA configuration 
corresponding to a given number of nucleotides iV per each strand. In that case, the rate 
equation for mutation becomes, 

2H 2 {A~ N ) = r (9) 

where r is the rate of mutation per nucleotide. In this case, the full degeneracy of the parent 
state as well as that of the daughter state have been taken into account. This leads to a 
very large magnitude of H even for a moderate value of N and reasonable value of r. Thus 
this violates the condition that H will have to be much less than A for the perturbation 
treatment to be valid. Therefore our formalism and analysis are appropriate for DNA se- 
quences corresponding to any definite species only. Moreover according to our prediction, 
the secondary repairing efficiency turns out to be of the same order of magnitude as the 
single mutation (point mutation) rate itself. This interesting result will have to be verified 
from the observational data after filtering out the tautomeric mutational component and 
taking into consideration the rate of production of the total 4N grand daughter configura- 
tion states generated after the 2nd round of replication from the singly mutated daughter 
states. This comparison will also throw some light on the validity of our modelling as well as 
the calculational scheme for studying the process of complementarity restoration. It may be 
remarked that amongst these 4N states, the one half of them would be the original parent 
state and the other half the base pair transformed states. Therefore we have totally (2N+1) 
distinct grand daughter states which are the products of the complementarity restoration 
process. 

4 Conclusions 

We have attempted a modelling of the replication accompanied by tautomeric mutation in 
DNA, by constructing a quantum spin-pseudospin based model analogous to some of the well 
known models used in statistical mechanics and magnetism [3-7,10-12]. The idea has been 
to incorporate quantum fluctuations to the existing ordered structure in a double-stranded 
DNA and calculate the rates of various processes. In particular, we considered both the 
disordering (mutation) as well as the order restoring CR events and calculated the efficiency 
of the restoring processes. It should also be highlighted that our result regarding the RE, 
is independent of the numerical magnitudes of the intrinsic coupling constants appearing in 
our model [9] . 

We have tried to apply our formalism and calculations to an experimental biological 
system. However, due to non-availability of sufficient data the comparison of observational 
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results with our prediction remains incomplete at present. To summarize, in our work we 
have chalked out a scheme for a possible novel repair mechanism for the DNA mutations, 
which is based entirely on quantum fluctuations. This is distinct from the conventional 
processes viz. primary repair system and the post-replication repair mechanism which are 
enzyme initiated [1] . It should be made clear however, that our scheme is applicable only to 
the case of tautomeric mutations. 

The possible future experiments should consist of detailed rate measurement of the pro- 
cesses of mutation as well as complementarity restoration arising purely out of tautomer- 
ization. For this the conventional methods like that of Luria and Delbruck [1,13] may have 
to be modified. This should enable one to critically examine our modelling and prediction 
and compare our theoretical results more precisely with the relevant observations. On the 
theoretical side too, there are lot of scopes for improvement and making the model more 
realistic by taking into account the intrinsic intra-strand couplings as well as the site and 
strand dependence of the transverse field. Last but not the least, this model suggested by 
us, although partially borrowed from condensed matter physics, has a lot of potentiality to 
describe the microscopic processes occurring within a DNA. This attempt is further inspired 
by our previous success with the analysis of the scaling behaviour of DNA involving the 
nucleotide distribution, by the methodologies of the statistical physics [14]. 
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